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Abstract 

The structure of polyhomogeneous space-times (i.e., space-times with metrics which 
admit an expansion in terms of log* r) constructed by a Bondi-Sachs type method 
is analysed. The occurrence of some log terms in an asymptotic expansion of the metric 
is related to the non-vanishing of the Weyl tensor at Scri. Various quantities of interest, 
including the Bondi mass loss formula, the peeling-off of the Riemann tensor and the 
Newman-Penrose constants of motion are re-examined in this context. 



1 Introduction 

In general relativity an important question is: what does the gravitational field of a radiating 
asymptotically Minkowskian system look like? The answer to that question proposed by 



Bondi et al. [|[, Sachs |29] and Penrose |27] seems to have been adopted by researchers (c/. 
e.g. ^), in spite of the wide evidence against this proposal: indeed it has been suggested 
both by the analysis of Christodoulou and Klainerman |^ and by various approximate 
calculations (c/. e.g. |jl2| and references therein) that such systems generically do not satisfy 
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the Bondi-Penrose-Sachs asymptotic conditions. In a recent study (c/. also j^]) of 
the asymptotic properties of solutions of constraint equations on spacelike hypersurfaces 
intersecting "Scri" transversally it has similarly been observed that generic Cauchy data 
constructed in such a setting by the "conformal method" failed to be smoothly extendable, 
after appropriate rescalings, to the conformal boundary. More precisely, it has been shown 
(c/. 1]^, ^ for more details) that, when considering Cauchy data constructed by the conformal 
method with smooth up to boundary "seed" fields and with the condition tr K = const 7^ 0, 
one has: 

1. generically, for such data no conformal factor exists for which the shear of X in the 
metric Q'^j vanishes at 5S; by the vanishing of the shear of I we mean the somewhat 
stronger statement that 

(1.1) 



X 



(recall that in the case of C^(A^) metrics the existence of an such that (|l . l| ) holds 
follows from the vacuum field equations, c/. e.g. [pl|]); 

consider those data for which the shear of X vanishes for appropriately chosen O. 
Generically, for such data the Weyl tensor of 17^7 dots not vanish at X. (Recall that 
for vacuum metrics 7 such that 0^7 is up to boundary on hA the vanishing of the 



Weyl tensor of 7 at X follows by a theorem of Penrose M, 28, 



The results obtained in seem to indicate very strongly that a consistent setup in which 
the gravitational radiation field can be described in generic situations is that of manifolds 
(A^,7), 7 = 17^7 with metrics 7 which are not smooth but polyhomogeneouM near X. (A 
function / is called polyhomogeneous if it admits an expansion in terms of log* r rather 
than , cf. Appendix^ for a more precise definition.) The object of this paper is to show 
that at least part of the results described above can be obtained in a rather simpler way in 
a Bondi-Sachs type setting, as set out in earlier papers in this series (Papers VII, VIII and 
IX, referred to here, frequently, as ||5|, ||2^, and Q respectively). 

In section 2, we show that the hypothesis of polyhomogeneity of X is formally consistent 
with the Einstein equations. (Note, however, that thanks to the important theorems of 
Priedrich |13, 14, 15 1, together with the results of [|l|, ^, a large class of space-times satisfying 



the Bondi-Penrose-Sachs conditions is now known to exist. On the other hand no proof 
that the Cauchy problem is well posed for polyhomogeneous but not smooth initial data of 
"hyperboloidal" type is available yet.) We show that the characteristic initial value problem 
of Bondi-Sachs type is formally well posed in the space of polyhomogeneous metrics, in the 
sense that the (retarded) time derivatives of the fields on the initial data hypersurface are 
polyhomogeneous if the free initial data are (with the same "degrees of polyhomogeneity" 
when these degrees are chosen appropriately, cf. Section || for details), and, in a manner 
completely analogous to that of the original Bondi-Sachs analysis, that one can write down 
a hierarchy of evolution equations for the coefficients of the polyhomogeneous expansion of 
the free data. 

In section 3 we show that in the class of space-times considered in this paper the con- 



formal factor can always be chosen so that (1.1) holds. Thus, Cauchy data incompatible 



^The term polyhomogeneous seems to have been adopted in the mathematical literature for the kind 
of expansion considered here, cf. e.g. [^. Gelfand and Shilov jl6| use the term associated homogeneous 
for a similar notion. The members of the Garchingrelativity seminar have suggested use of the term 
polyloganthmic for this kind of expansion. Winicour |32] uses the term logarithmic asymptotic flatness in a 
somewhat similar setting. 
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with (|1.1|) cannot lead to a space-time of the type considered here (c/. Q for a similar result 
in a somewhat different setting). We prove that initial data, constructed by a Bondi-Sachs 
procedure starting from free data smooth at T, will be smooth at X if and only if the free 
initial data are such that the Weyl tensor of 17^7 vanishes at X. We find that the Trautman- 
Bondi mass loss formula |o|, |, H] remains unchanged in the polyhomogeneous case; thus 



the 
timi 



ondi mass is still well defined, and is a monotonically decreasing function of retarded 
. We also note that for a class of polyhomogeneous metrics the peeling~off property of 
the Riemann tensor is the same as the one for smooth metrics up to 0(r~^~^), < e < 1, 
terms. We show that some quantities built out of the restriction of the Weyl tensor to X 
are (in general nontrivial) constants of motion, as already noted by Winicour |^2| and by 
Christodoulou and Klainerman Q. More generally, we find (in section 2) that the "leading 
log coefficients" of the polyhomogeneous expansion are constants of motion. We argue, 
from an explicit calculation in the axisymmetric case, that the Newman-Penrose constants 



of motion |23, 28| cease to be constants of motion in generic polyhomogeneous situations, 
although our example does give a new constant of the motion. (It could, therefore, be that 
some new functionals of the field, which reduce to the Newman-Penrose constants of motion 
when X is smooth, are constants of motion in the polyhomogeneous situation. We do not 
have an answer to that question.) Section 4 considers the construction of Bondi coordinates 
in our more general setting. 

The results of our analysis show that the presence of some log r terms in an asymptotic 
expansion of the metric is quite natural, and does not lead to any serious extra difficulties 
in the analysis of the geometry. Recall that the imposition of the conditions which lead to 
the vanishing of the log terms was interpreted in some earlier papers of this series as an 
outgoing radiation condition [^, 29, Two concerns had to be addressed: the possibility 



of advanced rather than retarded solutions, and the possibility of retarded waves travelling 
in the inward radial direction but at indefinitely large distances. With the help of our 
present understanding of Scri, it is clear that if Scri"*" is well-defined, as it is here, there is 
no advanced wave involved, and that a space-time has purely outgoing radiation if and only 



if there is no radiation at Scri~ ( cf. also 1 2C ] for a similar point of view and for explanation 
of the difficulties that arise with local characterization of incoming and outgoing parts of 
the field even for linear theories in fiat space). Since the space-times discussed here can 
satisfy both these requirements, we can safely abandon the "outgoing radiation" condition 
of §. 

Moreover, we note that there exists a family of electrovacuum "small data" space-times 



constructed by Cutler and Wald [11 1, and also a family of "small data" Einstein- Yang-Mills 
spherically symmetric space-times constructed by Bartnik Q, which have the following 
properties: they possess a smooth Scri''" and a smooth Scri^, and decay to a smooth i~^ in 
the future and a smooth i~ in the past. Because the metric decays smoothly both in the 
future and the past there is both outgoing and incoming radiation in those space-times. 
Since both Scri's are smooth, the "outgoing radiation condition" holds at Scri"*", and an 



^More precisely, for all polyhomogeneous metrics there is a quantity which we call the Bondi mass, which 
is a nonincreasing function of retarded time, and which reduces to the quantity defined by Bondi when 
Bondi's hypotheses are satisfied. We believe that the "real mass" should not be defined ad hoc, but by 
a limiting procedure involving perhaps the Freud "superpotential" for Einstein's energy, as done e.g. by 
Trautman in js^]. If one does that, we expect that one will find equality of the quantity we define as the 
Bondi mass with the quantity obtained from the limiting procedure only when V has no log^ terms, i.e., 
when the logarithmic terms in V start at the r~* level, with some i > 1. A precise formulation of such 
statements lies outside the scope of this paper. 



3 



analogous "incoming radiation condition" holds at Scri^, which is clearly absurd. These 
examples show that not only is the vanishing of the log terms at Scri unnecessary to ensure 
"outgoing radiation", it also does not prevent incoming radiation. We conclude that the 
association of the absence of logr terms with "outgoing radiation" lacks justification, and 
that a better formulation of an "outgoing radiation" condition could probably be given in 
terms of constancy of the Bondi mass at Scri~ . 

In the discussion above we have adopted what we consider to be now a standard notion of 
"incoming" and "outgoing" radiation. In particular, it is clear, from the hyperbolic nature 
of the Einstein equations, that whenever a conformal completion of the space-time exists 
in which the conformal boundary is an incoming null topological surface, then there can be 
no influx of gravitational radiation (or, for that matter, of any non-tachyonic matter fields) 
through the surface in question. Thus, the existence of a conformal completion of the above 
described nature guarantees that we have an isolated system evolving in a^elf-consistent 
way, regardless of whether or not the fields are asymptotically MinkowskianH, regardless of 
the decay rates of the fields towards the conformal boundary, the degrees of differentiability 
of some perhaps conformally rescaled fields at the conformal boundary, etc. In view of that 
observation it might not be so surprising that for the polyhomogeneous Scri's considered 
here the Trautman-Bondi mass-loss law holds, regardless of the occurrence of some perhaps 
high powers of logr in the 1/r terms in the metric. 

It should be pointed out that several of the results discussed here have already been 
observed in a similar setting by Winicour [^. (However, we learned about Ref. [^2| only 
after most of the work presented here was completed. Also it seems that in [32| emphasis 



is put on somewhat different issues.) In this context one should also mention the results 
of Novak and Goldberg (c/. also [10, |2^), who perform a somewhat similar analysis of 



the Newman-Penrose equations on a null initial hypersurface. 

2 The Bondi— Sachs characteristic initial value problem 

In this section we shall consider the initial value problem for space-times (7W , 7) with a 
metric of the form 

T/p2/3 

7^^ dx^ dx" = —du'^ - 2e^^ du dr + hab{dx'' + U''du){dx^ + U^du) . (2.1) 

We shall mainly be interested in the behaviour of 7 on the hypersurfacei 

= {u = 0, r > i?, G ^ 

where S"^ is topologically a two dimensional sphere. (The question of the existence of 
coordinate systems in which an asymptotically Minkowskian metric takes the form ( ^.ij) is 
considered in Section ^.) As has been analysed by Bondi et al. |5| in the axisymmetric case 
and by Sachs in general (c/. also P]), to construct a vacuum metric of the form (|2.1|) one 



^As pointed out below, several results proved in this paper (in particular the self-consistency of the 
polyhomogeneous setup) will still be true if the "sphere of null directions" is replaced by an arbitrary 
two-dimensional, perhaps but not necessarily compact, manifold M^. Examples of vacuum space-times 
with such an asymptotic structure (and actually a smooth Scri) are given by e.g. some Robinson-Trautman 
space-times. 

*Most of the analysis presented here goes through when 5*^ is replaced by any two dimensional, compact, 
orientable manifold ^M. 
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has to prescribe on N the family of metrics h{r) = hab{r, x"") dx"' dx^ on 5"^ parametrized 
by r, the family of vector fields U{r) = U'^{r, x"')da on S"^ parametrized by r, and the scalar 
fields V and (3. These quantities are not freely specifiable but have to satisfy constraint 
equations: 

V X G TAA R^^k^'X" = 0, (2.2) 

where is any null vector field tangent to M {e.g., k'^dfj, = dr). As has been emphasized 
in Q and [|^, the equations ( ^ ) do not impose any restrictions on hah dx°' dx^, and in fact 
can be viewed as equations which together with appropriate boundary conditions determine 
V , (3 and U"" da given hah dx"" dx^ . In |2^ it was shown that if we assume hah £ C°^{M) 
and moreover 

hAr. X") = HJ^) + + ^<^2^ + O (^) . (2.3) 

for some functions hab{x°-), h\f^{x°-) , a{x°-) , then we wih obtain r~'^V, (3, U°-, € C°°{M). 

In H, ^] the absence of trace-free terms in (2T^) was termed the "outgoing wave 
condition". This condition was imposed a priori in |5|, because the occurrence of the 
trace-free terms led to r~^ log* r terms in U"", V and subsequently in — this in turn led 
to r~^ log* r terms in hat at any later moment of time. It is therefore clear that a correct 
setup for analysing the characteristic initial value problem for metrics of the form ( |2.1| ) is 
that of metrics hah dx"" dx^ which are poly homogeneous to start with {i.e., admit an asymp- 
totic expansion in terms of r~^ log* r). It is the aim of this paper to re-examine both the 
constraint and the evolution equations for metrics of the form (|2.l] ) in a poly homogeneous 
setup. 

Before proceeding to a detailed analysis of the Einstein equations, let us first consider 
the question of the boundary conditions satisfied by the fields under consideration. Let 
therefore a polyhomogeneous metric hat dx'^ dx^ (see Appendix ^ for precise definitions) be 
given, and suppose moreover that hab S C*^(AA) (if we write hab G A^^'^^ , then the hypothesis 
hab S C^{M) is equivalent to the condition Nq = 0). It follows that the limits 

hah = hm hab 

r— >oo 

exist, with hab £ C°°(S^). As is well known, (c/. e.g. Q for a simple and elegant proof) there 
exists a diffeomorphism $ : S"^ ^ 5^ such that we have (t* h = cjP'h, where < (/> G C°°(5^) 
and h is the standard round metric on S*^, 

hab dx^ dx^ = de'^ + sin^ 9 dip'^ . (2.4) 



Replacing {r,x°') by (r,x°) = {cl}r,^°'{x^)) one obtains a metric of the form ( |2.1| ) in which 
(dropping bars on f, x") 

hm hah dx" dx^ = dO'^ + sin^ 6 dip^ . (2.5) 



It is not too difficult to show from eqs. ( |2.2D (which are written out in detail in Appendix 
^ that under the condition hah G Aphg the limits 



H = lim /3, (2.6) 
X" = - lim (2.7) 
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exist, with H, X"" G C°°(5^). Suppose for a moment that there actually exists a space-time 
with a metric of the form (2.1) on a set [7^ = {u € (— e, e), r > R, x"' £ S^} with some 



e > 0. (We should stress that in all the results obtained in this section the hypothesis of the 
existence of an evolution of the initial data defined on f/^ for some e > is not necessary. 
This is due to the fact that all our analysis involves only equations on M). Let ip{u,x"') be 
the one parameter family of diffeomorphisms of S"^ generated by the vector field X = X"" da] 
we thus have 

It is easily seen that the coordinate transformation {u,r,x^) [u^f^x"") = (n, r, x"), with 
implicitly defined by = 'ip'^{u,x^), leads to a metric of the form ( p.l| ) for which we 
have 

lim = . (2.8) 

r— >oo 



The above argument shows that the vector field X"" da defined by (2.7) has a gauge character, 
at least from a four dimensional point of view. It should also be pointed out that the 
transformations leading to (^]^) and (^^) are compatible with an initial value setup, because 
they do not deform the initial hypersurface M in space-time. This shows that there is no 
loss of generality in assuming that (|2.5| ) and (^) holdu. 

It is natural to ask whether H defined by ( p. 61 ) can be removed by an appropriate choice 
of gauge. As shown by Bondi et al. ^ in the smooth axisymmetric case, and as shown in 
Section Q in the general polyhomogeneous case, the condition 

lim P = (2.9) 

can always be achieved, at the price, however, of deforming J\f in space-time. Since it is our 
goal to analyse an initial value problem in which a null hypersurface M is given, we shall 
not assume that (^) holds unless explicitly specified otherwise.! (We will, however, find 
it useful to impose ( 2^ when discussing the physical properties of the four dimensional 
space-time. This is clearly justified by the results of Section ^.) 

Let us show that poly homogeneity of hah implies that of (3,U"',V and 

Proposition 2.1 .• Given any sequence {Ni}^^, Nq = 0, there exists a sequence {A^jj^g 
with No = 0, Ni = Ni, Ni > Ni, such that for all hab G ^^^'^ n C^{N) satisfying 
lim^^oo dx"" dx^ = dO^ + sin^ 9 dip"^ we have 

1. 

P, ^7^ r- V € Aphg n CO (AT) ; (2.10) 

2. If moreover 

lim U" = (2.11) 



^More precisely, there is no loss of generality in assuming that (2.5) holds, and there is no loss of generality 
in assuming that ( ^.^ ) holds provided the time-development of the data is sufficiently regular asymptotically 
in some neighbourhood of A/" (e.g., polyhomogeneous) to be able to perform the construction which leads to 



Cf. also [ p^ , for an analysis in which ( |2.9| ) is not assumed to hold, motivated by rather different 
considerations. The function here corresponds precisely to the function po = limr-,oo p of Il9|. 
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holds, then we have, for any j > 0, 

r-'v eAphgnC'>{/^), (2.12) 

|^y/i,fee^W>ncO(Ar). (2.13) 



Remarks. 1. It is rather clear from the gauge character of hm,._^oo that ( 2.1l| ) is not 



necessary for (2.13) to hold, with possibly a different sequence {Ni}. We have indeed verified 



this explicitly in the axisymmetric case (the validity of (|2l^) is in this case guaranteed by 



a cancellation of some "dangerous terms" which occur in the equation for ^i^)- It should 
be noted that, assuming the initial value problem for polyhomogeneous data on J\f can be 
solved, the u-dependent terms arising when limr_+oo f/" 7^ would lead to a u-dependent 
limit of hab differing from the round metric on surfaces other than M. This can also be 
seen from inspection of the transformation used above to remove X". It does not affect the 
Proposition above which applies only on M. 

2. Assuming the characteristic initial value problem for polyhomogeneous data can be 
solved in the space of space-times with a polyhomogeneous Scri, one way of understanding 
the significance of the relation between the {Ni}^Q and {Ni}^Q sequences is that if we 
are given initial data in a space characterized by the sequence {Ni}^Q, then the sequence 
{Ni}^Q can be chosen to be the one appropriate to the evolution of that initial data. 
Sequences {Ni}'^Q characterize those spaces which are invariant under evolution governed 
by the vacuum Einstein equations. 

Proof: Replacing hab by hat and r by r ^"^^^ g"'' we may without loss of generality 

assume \/ det h^t = sin 9. A simple but somewhat tedious analysis of the equations of 
Appendix making use of Proposition AJ in Appendix gives the following: the limits 



lim^^oo /5) lim^^oo U"", exist and are, respectively, a smooth function and a smooth vector 
field on a sphere. We define 



H= lim 13 , X" = - lim U" , 

r— »oo r— »oo 

ij" = 2e^^V''H , (2.14) 

where is the covariant derivative with respect to the metric h = limr^oo hab dx^ dx^ on 
5^. We then have 

P-He-^, (2.15) 

U^ + X^-^€^, (2.16) 



Let us also define 



r~' V G Aphg n C"(AA) . (2.17) 
ij = e^^ (1 + 2Af^H + 4|Pi/||) , (2.18) 



where is the Laplacian of the metric h, and | • |^ denotes the norm in h. If moreover 
lim^^oo U"" = 0, it follows that 

V-i/jre Aphg , (2.19) 
dhgb ^ C°°(5^) ^ Aphg 



+ (2.20) 

ou r 
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In the above analysis the only not entirely trivial step is to prove (|2.20|) . Indeed, equations 
(C.5) and ( |C.6| ) take the form 

"'- + /</'2 = Ci, (2.21) 
/</>! = C2, (2.22) 



dr 
d4>2 



dr 

where, using the notation of and Appendix j^, 

61 = rcosh(2(5) — , 
ou 

86 

(P2 = r — , 
ou 

f = 2sinh(2<5)^, 

and the Ca, a = 1, 2, can be found in Appendix T he hypothesis hab G n C°(7V') is 
equivalent to 7, (5 G A^^^ n C°(AA), so that from (pl5|)-(^7il|) one finds 

r/, Cae^, a = 1,2. (2.23) 

It is an exercise in ODE'S to show that under (|]2|) any solution of (|]2l|)-(|]2|) IS neces- 
sarily polyhomogeneous, with 

[The result can be proved by e.g. setting up a contraction principle argument in weighted 
spaces of the kind used in Appendix |^.] This implies G r^^ (7°° (5^) + r^"^ Aphg, as 
claimed. 

To prove our claim about existence of self-consistent sequences {Ni\ which have the 
property that hab ^ A^^^^ is formally preserved under time evolution, define 

N = Ni, 

so that 

\hab-Kb\=0{r-^ log^r). 

Consider a term, say Xi in hab which has a radial behaviour log-* r. From the equations of 
Appendix P one easily finds that such a term produces terms r~*~^ log^"'"-' r + lower ordeiQ 
in /3. Next, if i = 2, such a term will produce terms r~^ log-^"*"^ r + lower order in U""., while 
for z 7^ 2 it will lead to terms r^*^^ log-' r + lower order in U"". Using the Kronecker 5\ 
defined as usual by 

, ^ c-h f 1, for a = 5, 

a,b ^ ]R , 



0, otherwise, 

we conclude that x produces terms r^*^^ log-'^^* r + lower order in U"". There is a cancella- 
tion in the equation for V which implies that if z = 1 and A^(= j) = 0, x will generate terms 
r^^ log^^ r -|- lower order in V , rather than the log r -|- lower order terms which would have 



''By "lower order" we mean terms which have the same power of r ^ and smaUer powers of logr, or 
higher powers of r^^. 
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appeared if the cancellation had not taken place: here A^i is determined by N2 (Ni = 
if = 0). In general, the cancellation implies that the leading term in V arising from 
i = 1 and A*" > is of order log^ r, while the leading contribution^ from x is of order 
log-'"'"''' r. Inserting all this information in the equations for one in general obtains 

a contribution r~^~* log-'"'"''' r + lower order from but an i = 1, iV > term contributes! 
only a term of order log^^^ r. 

To construct a sequence {Ni} given a sequence {Ni}, set A^i = Ni. Then from the 

du 



equations for one obtains, if A^i > 0, 




O (r-^ log^i-i r) , (2.24) 



for some h^i){u, x"'). If A'^i = 0, the next contributions to the u-derivatives of the 7 and S used 
in Appendix 1^ are 0{r~^ log^^"*"^ r), which implies that the time derivative of the trace-free 
part of hab — h\i^{u, x")/r is 0{r~^ log^^"*"^ r); on the other hand, the time derivative of the 
"trace part" of hab — h\ij{u, x°-)/r is 0(r~^) and is determined uniquely by the ^-derivatives 
of h\j^{u, x°). ( p. 24 ) shows that the coefficients of the log* r, i = 1, . . . , A''i are constants 



of motion, and if we set N2 = max(A''2,A^i — 1) then the space A^^"-^ will be formally 
preserved by evolution up to 0{r^^^''), e > 0, terms. Proceeding recursively one can 
construct a self-consistent sequence {Ni}. The analysis of the higher n-derivatives proceeds 
in a similar manner by considering the equations obtained from the Einstein equations by 
M-differentiation, and the result follows. □ 
From what has been said in the proof above it should be clear that if we write, along 

00 Ni 



hab '^^^hijabiO,ip)r 'log 



r 



i=0 j=0 

then on M from the Einstein equations one obtains a Bondi - van der Burg - Metzner type 
hierarchy of equations 

Pi jab 



dhijab 



du 

where Fijab is a function 0, if and the hkicd-, < A; < i — 1 together with a finite number 
of their derivatives. In particular, if we assume that a poly homogeneous expansion of the 
metric also holds in a neighbourhood of N one obtains: 



Proposition 2.2 .• Under the hypotheses of Proposition \2J\, including (2.11), we have: 



1. The coefficients of r ^ log-' r, 1 < j < Ni, in hab f^i"^ constants of motion. 

2. Suppose that Ni = 0. Then the coefficients of r~'^ log-' r, 1 < j < N2 = N2 in hab CLnd 
the coefficients of r^^ in the trace-free part of the hab ore constants of motion. 



*The radial behaviour ' log-'^** r is obtained here by a sheep calculation. It seems that knowledge of 
some cancellations for i = 1 is necessary for the argument to hold. Nevertheless for i > 1 a straightforward 
analysis of the V equation (C^) yields a leading order contribution r^^* log-'^''* '^^i r to V from without 



going into the details of the cancellation structure of the equations (which sheep automatically does). One 
cou ld use this estimate of the contributions of x to 1/ in the remainder of the argument to prove Proposition 
with a perhaps somewhat "worse" sequence Ni. 
^This exponent is, again , obtained using sheep. Analytically it is more or less straightforward to estimate 



the right hand side of ( 2.24 ) as 0(r ^ log^^ r ) . Such an estimate would be sufficient for the main conclusion 
to remain valid. 
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3. Suppose that on M we have 

hl.ix'') , a{x'^)hab{x'') 



h., - lim h^, = + ^H^iWi ^ o (.-3 log-3 r) , (2.25) 



with hab given by ^2.^ , and that Ni = . . . = Ni = 0, i > 2. Then we can set Nj = Nj, 
j = 0, . . . ,i + 1 and the coefficients of r"*"-*^ log-' r, 1 < j < Ni^i = Ni+i in hab «re 
constants of motion. 



Remarks: In the notation of |^], ( 2.25| ) is equivalent to 7 = c/r + 0{r ^log ^r), 



5 = d/r + 0(r~3 log^^ r). Note that the time dependence of the contributions to hab is 
undetermined by the above considerations, which are purely asymptotic0. In Section 4 we 
show that the only remaining coordinate freedom is given by the BMS group, which implies 
that we have only one arbitrary function of 9 and (j) available for changing the values of the 
constants of motion just obtained, so that not more than one of them (if any) can be set to 
a fixed value. 

It is natural to consider extending the second result in the Proposition and ask what 
is the "smallest" self-consistent sequence Ni if A^j = for all i and it is not assumed a 
prion that (|]2|) holds. When hm^^oo = 0, one easily obtains, from what has been 
said earlier, that we can set 

iVo = iVi = iVa = 

and 

iV, = 1 . 



Moreover in such a case the arguments of the proof of Proposition 2J show that we will 
have 

fi e ^^^'>, n^ = n^ = n^ = nI = q, ivf = i, 

V-r e ^i^n, N^ = Q, nY = 1. 
In asymptotically Minkowskian coordinates 

(t, X, y, z) = {u + r sin cos r sin sin (/), r cos 6) 
this corresponds to a metric ( |2.lD which along M approaches the Minkowski one as 



7^. - = ^ + ^^^^ + ^ + 0(r-'-^) , e>0. 



The results of Proposition |2.1| can be generalized to include both matter fields and 
rather weaker asymptotic conditions. Let us start by defining a space of functions C^^: for 
/i, A G iR a function / will be said to be in C^'^(AA), or for short in Cf^^, if for all i G W 
and for all multi-indices a we have 

|«/| <a,ar-''-*(l + |logr|)^ , 



^°These contributions define an analogue of what is called a "news function" in and should be 

determined by the behaviour of sources and/or of the gravitational field in the interior, or perhaps by some 
interior boundary conditions, in a complete solution. 
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for some constants Cj^c, where the coordinates v stand for 6, <j). In order to avoid in what 
follows a rather annoying discussion of some not so interesting special cases we shall always 
assume that in all spaces considered if // = 1, 2, 3 or 4 then the logarithmic behaviour 
exponent A satisfies A 7^ —2, —1. We wish to show that the space of metrics on J\f of the 
form 

hab - Kb G r-^AP^3 + C7^\ > (2.26) 

is formally preserved by evolution with the Einstein equations with possibly some matter 
fields. Let us denote by {l.h.s.}^, respectively by {r.h.s.}^, the left hand side, respectively 
the right hand side, of the n'th equation of Appendix |^. Then, in the presence of matter 
the Einstein equations become 



{l.h.s.li = {r.h.s.li + -Tn , 

{l.h.s.}2 = {r.h.s.jg + 2r^Ti2 cosec 9 , 

{l.h.s.jg = {r.h.s.jg + 2r2fi3 , 

dV e^'^/i"'' - 
— = {l.h.s.}4 = {r.h.s.}4 — Tab , 

{l.h.s.}^ = {r.h.s.ls + — (e-2^r22 - e2^sin-2 0T33) 



{l.h.s.}6 = {r.h.s.}6 + 



2r cosh 25 sin 6 



'23 



-K'^Tcdh23 



where h""^ is the matrix inverse to h, 



(2.27) 

(2.28) 
(2.29) 

(2.30) 
(2.31) 
(2.32) 



T 



K is the gravitational coupling constant, and T^^ is the energy momentum tensor of the 
matter fields. We shall require 



fia G r-2^P^9 + C7^i-'^i° , a = 2, 3 

~2^T22 G r-^AP^3 + C^r,A. 
^23 - \K''f,dh23 G r'^AP^B + C^.-A. 



sin~2 6»f33 



^11 > 2 , 
/U12 = /il3 > 1 , A12 
/"o > , 



\13 



(2.33) 
(2.34) 
(2.35) 
(2.36) 

(2.37) 



where the various powers in front of A^^^ and the exponents in C^^''^' have been chosen 
so that the leading order behaviour of the various functions which appear in the metric 
coincides with the behaviour one observes in the vacuum case. Assuming that limj.^00 P = 



lim. 



r->oo U"" = 0, an analysis as described in the proof of Proposition 2.1 leads to 



■^Kb G r-^p^s + c; 

ou 



00 + f-^oo 



It follows that for 



H < min{^ii - 1, /xia, + 1, /io + 2} 



(2.38) 
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(with appropriate inequalities for the A=k 's if the above inequahty is an equahty) the Einstein 
equations will formally preserve the space of metrics hab satisfying ( 2.26| ). The proof of ( 2.3S| ) 
follows immediately by integration of the equations ( |2.27D ~( ^.32 ), and of the equations which 
are obtained by ^-differentiation of those. Step by step one obtains: 



V 



■a 



„/Jla — 2,Ala 



Mia 



[If we 



the above equations inserted in the evolution equations ( |2.31 )-( 2.32| ) yield ( |2.38|) 
assume that the appropriate decay properties are also satisfied by d^Tf^^/du^ on M, j = 
0, . . . , J, then a preliminary examination suggests that d^hab/du^ will also be of the form 
(|2]3|) for i = 0, . . . , J.] We note that the previous calculations on similar lines p6| 



led to upper bounds on the parameter /i which do not appear here because, unlike those 
previous calculations, we do not forbid the appearance of log r terms. 



3 Geometric Interpretation 



Consider a metric of the form ( |2.1[) . Following Penrose |27] it is useful to introduce a new 
coordinate 

X = , 

so that by Proposition ^ when hab G -A^^^ n C°(A/') the metric 

7^,y dx'^ dx'^ = x^ jfj^u dx^ dx^ 

= -Vx^e'^'^du^ + 2e'^^dudx + hab{dx^ + U^du){dx^ + U^du) (3.1) 

is polyhomogeneous on the set {x G [0, l/i?], G S"^}, i.e., there exists a sequence {iVj} 
and functions ^fiuij 

oo TVi 
i=0 j=0 

No = 0. 

(When (|2.25| ) holds and hab £ C°°(A/') one actually obtains iVj = for all i, and the metric 
( |3.1D is in C°°(Ar). This corresponds to the standard Bondi-Penrose-Sachs situation of a 
smooth Scri.) We have the following result, which is established by calculating the Christoffel 
symbols of the metric x'^j^i, using sheep, and making use of vacuum Einstein equations for 
7^j/ in the first leading orders: 

Proposition 3.1 .• Consider a characteristic initial data set with hab £ -^hg ^ C^{J\f), 
limr^oo hab = dO^ + sin^ 9 dip"^ and lim^-^oo U"" = 0, and set = x = r^^. We have 

lim VaVfef] = , a, 6 = 2, 3 , (3.2) 

X—^0 

where V is the covariant derivative of the metric j^^, = x"^^^^. If moreover lim^^^o = 0, 
then we also have 

lim V„Vi.J] = , /x,z/ = 0, ...,3. (3.3) 
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Recall that the geometric meaning of ( |3.2| ) is the vanishing of the shear of the hypersurface 
{x = 0} (c/. e.g. [^]). Another interpretation of ( p.2p is that the conformal extrinsic 
curvature of Scri vanishes, cf. [Appendix] . This property of Scri is well known for 7^;/ G 
C°°{Af). It should be stressed that for general polyhomogeneous metrics as constructed 
(on Af) in the previous section, the fact that the left hand sides of (|3.2| )-( ^^ exist and are 
bounded is a non-trivial statement which makes use of vacuum Einstein equations to the 
first two leading orders, because equations (|3.2| )~( |3!3|) contain derivatives of the metric which 
could potentially blow up as log^^ x as x ^ 0. Proposition 3.1 is the original observation 
which led to the proof in |^], that generic Cauchy data constructed by the conformal method 
as in will lead to space-times which cannot admit a polyhomogeneous Scri. 



Throughout the remainder of this section we shall assume that 
hm habdx^dx^ = dO"^ + sm^ 6 dif^ =: habdx^'dx'' 



lim U" = 

r— >oo 



(3.4) 
(3.5) 



lim (3 = 0, 

r— ►00 

A textbook property of smooth Scri's is that the Weyl tensor of the conformally rescaled 
metric vanishes at Scri, cf. e.g. |31|. A sheep calculation of the Weyl tensor of the metric 
7^jy dx^ dx'^ , assuming that 7^1, is vacuum, gives: 

Proposition 3.2 .• In addition to the hypotheses of Proposition \3. j| , let hat G C'^{M), and 



let (3.4)--(3.5) hold. In local coordinates define 



Xab 



lim 



ab 



dx'^ 



2^ 9x2 



ab 



Let Cap-fS denote the components of the Weyl tensor of 7 in the half-null tetrad 

.1 x^V 



-du — dx , 



= — (cosh Se' ^ sin Q sinh (5 e '^U'^^ du + cosh S e'^ dO + sin 9 sinh 5 e ^ d(p , 
9^ = - (sinh 5e^U^ + sin 6 cosh 5 e'^U'^') du + sinh 5 e'<d9 + sin 6 cosh 5 e'^ dip . (3.6) 
Then we have 

lim Ciaib = Xab , (3.7) 
I— >o 

while all the remaining components of lima;^o Cap^yS vanish, except of course those which 
can he obtained by appropriate permutations of Ciaib- 



[Let us point out that Proposition 3^ will still be true with non- vanishing Ricci curvature 
provided that (|23^ ) and (|2.33h ([237|) hold with ^o,/^r > 1, fiia > 2, fi, flu > 3, and that 
the powers of r^^ in front of the polyhomogeneous pieces in T^^, in eqs. ( 2.33| )- (|2.37| ) are 
increased by one.] Proposition shows that the Weyl tensor of 7^1, vanishes at x = if 
and only if the trace free part of q^^'' \x=o vanishes, i.e., if and only if what Bondi et al. 
termed the "outgoing radiation condition", ( p. 25 ), holds. Let us also mention that ( ^.7]) is 
equivalent to 



x=0 



^-1 



x=0 



^5 



x=0 



x=0 



x=0 
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where the "S/j's are the Newman-Penrose components of the Weyl tensor of the metric 7 in 
a null tetrad related to the tetrad (p. 61) above in the obvious way. 



It is worthwhile emphasizing that it follows from Proposition 2^, point 2, and from (| 
that the components Ciaib(u , , x'^) are pointwise constants of motion, i.e., independent of 
u. This result seems to have been already observed by Winicour [^2|, and independently 
by Christodoulou and Klainerman [^] (under much weaker asymptotic conditions). 

Perhaps the most important result of the Bondi-Sachs analysis is the well known theorem 
(originally due to Trautman |3^) that Bondi's mass is a decreasing function of u. For the 
metrics under consideration here, let us define the Bondi mass as (—1/2) the integral over 
the sphere 5^ of the coefficient in the expansion of V. This definition clearly reduces to 
the original one by Bondi-Sachs, when the conditions imposed by Bondi and Sachs hold. We 
have found that, under (|3.4)~ (|3.5| ), and whatever the sequence {Ni} and the hab G A^^'^ , 
the mass loss equation (35) of which can be obtained by equating to zero the integral 
over of the right side of (|3.8|) , (c/. also [^], eq. (4.16)) remains unchanged. This can be 
seen as follows: under the above conditions it follows from the vacuum Einstein equations 
that 

V -r, r2/3, r'^U'' G AP^3 . 



As has been shown by Sachs |29| following the original observation of Bondi et al. in 



the coordinate system of (2.1) we have (if the other field equations hold) 

^00 = lim r'^Roo = . 

A SHEEP calculation gives 

hm r^Roo = lim -r^h'^'h"'' ^ ^ ^ ^ — , (3.8) 

r^oo r^oo \ (ju 4 ouor ouor ou 



where Da is the covariant derivative of the metric hab- It is clear from (3^) and from what 
has been said before that those log terms which could contribute to this equation, if any, 
drop out because their u derivatives vanish. It follows that for all polyhomogeneous initial 
data the Bondi mass is a non-increasing function of time whenlll limj.^00 /? = 0. 

We would like to point out that it is clear that the mass as defined above is the "correct 
mass" for polyhomogeneous initial data hab G A^^^^ with Nq = Ni = 0. On the other 
hand we believe that some care should be taken when interpreting the above as the mass 
in the case A'^i 7^ 0, since in that case the leading order behaviour of y — r is logarithmic, 
which might refiect an infinite or ill defined mass of the system. Such a possibility is also 
suggested by eq. ( |2.30| ) which shows that some log terms might arise from 1/r terms in 
h°-''fab- Now in an orthonormal tetrad in which e^ is timelike one finds that h'^'^fab = 
—Kr'^{TQ + Ti) = Kr'^{Too — Tn), so that 1/r terms in h°'^Tab correspond to an infinite 
amount of matter energy: if r^Tgo behaves as 1/r, then Tqq (the matter energy density) 
ceases to be integrable over M. A thorough analysis of the conditions under which the mass 
at null infinity is finite and well defined lies outside the scope of this paper. 

It is natural to ask what happens with the "peeling-ofF' property for space-times for 
which ( |2.3| ) fails to hold. A sheep calculation shows that (in the vacuum case) for any 
hab G C°(AA) n AP''9 along Af we have 

p _ ^al3"/S ^al3-y5 ^a^-fS ^ ^ajS-yS 



When Unirnex) /3 7^ 0, then the "mass aspect" is no longer given by the standard Bondi-Sachs expression. 
But monotonicity of an appropriately defined total Bondi mass of course remains valid. 
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with R^f^^g and R^ibjS peeling off exactly in the same way as they would if (|2. 



were 



satisfied, i.e. R^/^^s type N, and R'^f^^g is of type III; in fact i?^ 
exactly the same as in the case considered hy Sachs (c/. also [^0|l). This follows from 



a/375 and Rf^f^^s are 



the fact that all terms which contribute to R 



1 



and R 



are ^-differentiated, and those 



log terms which could potentially contribute at this order are constants of motion. Here 
R^ap-yS ™ay contain powers of logr. 



In the case 7^1/ G C°°{jQ') it was observed in [23, 24] that there exist some nontrivial 
global constants of motion for a vacuum gravitating system. We wish to point out that 
these quantities cease to be constants of motion even in the case in which hab £ C°°(AA) if 
one does not assume that ( |2.25| ) holds. Clearly it is sufficient to prove that assertion for 
those metrics (|2.lD which are of the Bondi-van der Burg-Metzner form |@]: 



, 



hab dx" dx^ = e^^ dO'^ + e'^^ sin^ 9 dip'^ , 
C/'^ = 0. 



(3.9) 

(3.10) 
(3.11) 



Let us expand the function 7 appearing in (|3.10[) as 



c , 72 , 73,1 logr 73 74,1 logr D 
7 = - + — H 3 h — H 1 ^ — + 

ry^ fj-i^ f^f-J ^f-J ry^^ ly*^ 



[The terms 73 1 and 74 1 above are necessary: even if 73 1 = 74 1 =0 (which we are 

free to assume), we shall have 73,1 / / 74,1 at later times in general, as a consequence 
of the evolution equations.] In the axisymmetric case with 72 = 73,1 = the constant of 
motion is given by [^ 



dfio 



sin dO ( 



If one takes the minimal sequence of Ni as described in section 2 and allows 72 ^ 0, one 
finds that the 73.1 and 74^1 terms (but not r"^ log-' r or log-' r with j > 1) arise. After 
a long calculation it turns out that the u-derivative of the 74^1 term, when multiplied by 
sin'^ 9, is a total derivative with respect to 9 and, removing further total derivative terms, 
one finds that the equation of motion for D takes the form 

^ = / Gd9 d(t) 
ou Js^ 



where 
G 



72 sm 
+15sin3 6l7^i |>/12, 



2q( a , 1^,^ in„. , o..:„fl„„.Z) 16 COS^ 6* C 



sin^ 6'(-4— 7 + 16M - lOu + 4c) - 24 sin 9 cos 9 

o9^ o9 



(3.12) 



Til 



73,1 



u=0 



and M is the usual Bondi mass aspect {i.e., —2M is the integration constant 



which appears when integrating the V equation (|C.4D ). 



15 



u=0 



For given M ^0, 73^1 (perhaps, but not necessarily, being zero) and c 



u=0 



seems obvious from ( 3.12D that the function 72 



u=0 



it 



dV 

du 



u=0 



u=0 



can be chosen so that we have 



It must be pointed out that the above argument falls short of being a rigorous proof: 
the function M has a global character, and in particular it might not be independent of 



72 



u=0 



and 73^ 



u=0 



(Nevertheless the above calculation shows that no obvious miraculous 



cancellations occur, and we find it completely implausible that in, say, the vacuum case 



there exists some kind of conspiracy between the functions 72 ; 73 1 , c and M which 

u=0 ' u=0 

leads to the identical vanishing of the ^-integral of G.) 

It is curious, and not entirely unexpected, that in the axisymmetric poly homogeneous 
setting there is a Newman-Penrose type quantity which is again a constant of motion. 
Define 



s 



£•2 



74^1 sin 6 d^Q . 



(3.13) 



As we show in Appendix Q is conserved by the evolution via the vacuum Einstein 
equations. It seems clear to us that an analogous result will be true in the general case, 
without assuming axisymmetry. 



4 Existence of Bondi coordinates 

Consider a metric 7 defined on the set 

^rAA = {r>R,ue ((7i, (72), r G (4.1) 

and suppose that there exists < e < 1 such that 

for {x^'x") = (iin),(nx"),(fx") 7£Mi- G Aphg, \jxf^x-\ < , (4.2) 

7™ e ■Aphg, e<lfu< e""^ , (4.3) 

7fr G r~'^Aphg, |r^7ff I < , (4.4) 

7£.i>fa G r'^Aphg, \r~'^lx''x-\ < • (4-5) 
Following Penrose |27] let x := f~^, and set 

7/J!^ — 7/^1^ ■ 

From ( |4.1| )-( ^^ it follows that the metric ^^^^^dx^dx'^ can be extended by continuity to a 
poly homogeneous metric on the set 

^i/rA,C2 = {0 < ^ < V^, ^ e (Ci, (72), x'^ G S'} . (4.6) 

Actually the conditions ([4.3|)-( |45| ) guarantee only that the appropriately rescaled functions 
7^jy can be extended by continuity to the boundary, with the appropriately rescaled metric 
degenerating perhaps at the boundary. We thus add the supplementary restriction that 



16 



"fabd'X'^dx^ is non-degenerate up-to-boundary, and that 7 is also non- degenerate with signa- 
ture ( — h ++) up~to-boundary, in a sense which should be clear from what is said below. 
Define 

X = {x = o,ue ((71,(72), r G S^} . 

Throughout this paper we shall suppose that T is a null hypersurface; as in the smooth 
case (c/. e.g. in the poly homogeneous case this will necessarily hold if 7 is vacuum. 

Let u be any smooth function on Z, and extend n to a smooth function defined in some 
neighbourhood of I in any way. Let w^dx^ be any smooth nowhere vanishing one-form 
field defined in a neighbourhood of Z such that Wf^X^ = for all E TI; from the fact 
that 2 is null it is easily seen that 



On X consider the one-form field 



dx^ = aw. 







dx^^ + du 



(4.7) 



(4.8) 



with some function a; from (|4.7|) it follows that the equation 



du 



) will have a (unique) smooth solution a\x provided that 



is bounded away from zero. 



(4.9) 



(4.10) 



which we shall assume to hold. (Note that ( 4.10| ) implies that = '^^^ky will be trans- 
verse to I.) By Proposition B.l point 1 for every ('u,x") G I there exists a null geodesic 
x^(s,M,i;") such that dx^/ds{0,il,x"') = k^\j . There also exists a diffeomorphism — > 
x"(x^) such that in the coordinates (u, x") we have at ii = 



lab 



u=x=0 



dx^dx^ = (j)'^habdx''dx^ = (p'^{de^ + sm'^9dip'^) 



where (p £ (7°° (5^) is uniformly bounded away from zero. Let (u, s,x°) be obtained by Lie 
dragging {u,x"') along the integral curves x'^{s,u,x"'). Thus if we set /c'^ = dx^/ds, then 



k''u,f,= 0, 



k^x ,fj_ — 0, 



n(0, n, x") = u , 



(4.11) 
(4.12) 



By point 2 of Proposition [B.l| and by the implicit function theorem there exists a neighbour- 
hood of T on which (n, s,x") form a coordinate system. As k'^ is tangent to null geodesies 
we have ^fj.uk^k'^ = 0, so that 

,d d , 

It follows that in these coordinates 7 takes the form 



^^lydx^dx" = ^uudu^ + 2%sduds + ^ab{dx°- + U''du){dx^ + U^du) . 



(4.13) 



17 



Setting f = 1/s the original metric 7 takes the form 

jf^udx^'dx" = r^^uudu^ - 2(^f%sdudf + r'^%b{dx'' + U"-du){dx^ + U^du) , 
where all functions in are polyhomogeneous, after appropriate rescalings. Finally define 



sin2 9 det 7^^ 



dx"" dx^ 
^x^^ dx^ 



(4.14) 



One easily finds 



dr log^f, 
— = (j) + 0[ — 
or r 



for some N , so that there exist constants R, Ci, C2 such that (r, u, x") as constructed above 
form a coordinate system on Ur^Ci ,C2 > where Ur^^ is given by ( 4J ) (in the coordinates 
{r,u,x"')). Going to this coordinate system from ( 4.14| ) one concludes that 

det(7ab) = r'^sin^6' . 



In this way one obtains a metric of the form ( |2.1D , which satisfies appropriate Bondi require- 
ments at n = 0. If one moreover assumes that 7 is vacuum, then the Einstein equations 
imply 

— ( lim %b , 



, 



so that we have 



lim 7ab = dO"^ + s\v? Od(p^ 

r— >oo 



for all u S (Ci, C2). As discussed in Section |2| we can achieve 





lim U" 

r— >oo 



by appropriately propagating the coordinates away from the surface li = 0. 

The above construction shows that the Bondi coordinates above are uniquely determined 
by the choice of a function u = u]^ . Let us show that the freedom in the choice of u| j can 
be considerably reduced if we require 



lim /? = 0. 



(4.15) 



To achieve ([4.15| ), let n|j be given by 

emu',x-)^^> + a(x'*) , {H = lim /?) , 



U[U, X 



(4.16) 



where a G C°°(5'^) is an arbitrary function. We can now repeat the construction described 
above of a coordinate system (s, -UjX") based upon this function u\j , obtaining again a 
metric of the form ( 4.13| ) in this coordinate system. Let thus a vacuum metric 7 of the form 
(2.1) be given, and set 7^1, = r~'^^^^, x = r~^, k^^ = j^^'^u,^', from the equality 



r 
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by construction of we obtain 



with 



-dx 

du 
dx 



-2H 



du 
dx 



du + Vu 



da , 



_ 1 

x~ 2 



1^-11 



where T> denotes the covariant derivative of the metric h on S"^ . We also have 



dx°' 



du 



dx"" 



dx^ 



Ob > 



dx 



V'u 



(4.17) 
(4.18) 

(4.19) 



Prom the "barred" equivalent of (4.14) 



sin^ 9 det 7^^ 



and from what has been said one easily finds 



dx"- dx^ 
dxt^ dx'' 



dr dr dr 

lim — = 1 , lim — = lim — — = 

r— >oo Qf r-^oo qu r^oo ox"' 



(4.20) 



Equations ( [4.18| )-( ^.2ClD show that in the coordinate system (n, r, x") the metric takes the 
desired form with 

lim /3 = . (4.21) 
It is worthwhile mentioning that the only freedom in the choice of coordinates left at this 



stage is that of the function a in ( 4.16 ). Any two coordinate systems {ui,ri,x1) and 
(■"2, ^^2) a^f) which satisfy our requirements will have the property that 



a = {ui- U2) 



is a function of x" only, and the coordinate system (M2,r2,X2) is defined uniquely by a 
and by (ni,ri,x^). It follows that for polyhomogeneous metrics the asymptotic symmetry 
group is the BMS group, as in the smooth case. 
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Appendix A Conventions, function spaces 

We shall assume that all manifolds we discuss (which will have dimension 2, 3 or 4) are para- 
compact, connected, Hausdorff and smooth. The summation convention is used throughout 
this paper. Space-time is as usual taken to be a Lorentzian 4-dimensional manifold with 
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signature (— , +, +, +) and metric connection r^^^, where Greek indices run from to 3; 
its Riemann curvature tensor is defined by 

and the Ricci tensor and scalar by 

Rgs = R^'paS, R = g^^R^u- 
With these choices the Einstein equations take the form 

with a positive constant k. 

In this paper we are considering asymptotic behaviour near null infinity. We will use 
M to denote a manifold with boundary so that M = int M is a manifold of dimension n 
and dM = dM is a manifold which will be assumed to have a finite number of connected 
components 9Mj. By an abuse of terminology M will also be said to be a manifold with 
boundary. As usual TpM will denote the tangent space to M at p; for p in dM one has the 
notion of "half-tangent space at p" which is defined in a natural obvious way, and we shall 
still write TpM for this space. 

Throughout the paper x will denote a defining function for 9M, i.e. a function satisfying 
^\dM = 0, X > 0, dx{p) / for p G (9M, and the implication x{p) = p G dM holds. 

We can always choose a finite number of coordinate charts (f)j : Oj = {y e : 

> 0}, j = 1, . . . , J, covering a neighbourhood of dM such that y^ = x. When referring 
to local coordinates we shall implicitly assume that y^ = x, and we shall use the letter v to 
denote the coordinates y^, . . . , y"-; 

i;^ = y^, A = 2,...,n. 

Thus y = {x,v) . The standard Schwarz multi-index notation is used throughout; thus if 
a = {ai, . . . an), then 

U — Uy — Uyn — — Uy , 

where (i = {02, ■ ■ ■ ,an)- 

For k G 1Nq° = M U {0} U {00} the spaces C]^p(M) are the spaces of functions fc-times 
differentiable on M. We have added the subscript "loc" to emphasize the fact that a function 
in Cjqj,(M) need not extend to the boundary of M (in this respect the subscript "loc" does 
not imply the same sense of "local" as the local coordinates we have just defined); similarly 
for A; > 1 even if the function itself extends by continuity to dM then its derivatives do 
not have to extend, etc. We use the symbol C''{M) for the Banach spaces of functions 
differentiable fe-times on M such that / and its derivatives up to order k can be extended 
to continuous functions on M, and equipped with the supremum norm. 

Let fi be a sequence of functions in C^^{M) and let xi > 0; we suppose that, given 
N G -2V, there is a sequence n — >■ 00 and some constants C, jv such that for all \a\ < N 
and for all < x < xi, 

\d^fi\ < Ci,Nx'''^ , 
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To express the notion of successive approximations good to all powers of x and for all 
derivatives in a precise sense, we shall write 



if for every n,m G IN there exists N G IN and a constant C{n, m) such that for all |a| < m 
and for < x < xi 

dy(f-f2f^) <C{n,m)x^. 
^ 1=0 ' 

Consider a sequence {A'jIj^q, Nj G Wq. / will be said to be poly homogeneous if 
/ G C^JJvl) and there exists a sequence of functions /jfc G C°°(M) such that 

oo A^j 
j=0 fc=0 

We write / G A^^^^ , and define ^p/^g = U{Af^}^^^^>. 

Remark: As formulated here, the fij may depend upon x. To avoid this would require 
the introduction of local coordinates near the boundary, a specialization which we do not 
yet wish to make. However, once we do fix a coordinate system, then we can Taylor expand 
each fij with respect to x to any finite order, so each fij has a poly homogeneous expansion 
with no log terms, and obtain an expansion for / with some (other) functions fij which 
depend only upon the coordinates v. 

A function f{r,v) defined on an open set of the form O = {{r,v) : r G (ro,cxD),i; G Q} 
for some suitable set Q will be said to be in C^{0), < < oo, if f{l/x,v) G C^{U), where 
U = {{x,v) : X G (0, l/ro),v G Q}. Similarly / will be said to be poly homogeneous on O if 
f{l/x,v) is polyhomogeneous near x = onU. 

Let -F be a function space over M. A tensor field X = (X'^p), where a, (3 are some multi- 
indices, |a| = r, \(5\ = s, will be said to belong to F if in local coordinates as described at 
the beginning of this section the components X°'i3 of X are in F. 

Let F be a function space. We shall write that f G r°' \og^ r F if r~° log~'^ r f G F. 

Let Fi, F2 be function spaces. We shall write that / G + F2 if there exist fa & Fa, 
a = 1,2, such that / = /i + /2- 

The following observations are useful when proving Proposition 



Proposition A.l .• 

fGr'^log'^rAphg, g e r'' log"" r Aphg =^ f 9 r''^'' log^-"" r Aphg ; (A.2) 

FeC^iM), feC^{Ar) + ^ =^ F{f) e A^hg ; (A.3) 

ieZ, fer' Aphg =^ Vi G W, r'-^ Aphg , (A.4) 

Va, d': f er'Aphg-, (A.5) 

ieZ, fer' Aphg =^ r Aphg . (A.6) 

JR 

Proof: ( [A.2| ) is easily proved starting with the following elementary observations: f,gG 
Aphg =^ f + g £ Aphg] f G log'^ r C°°{dM) (note the Remark above), g G Aphg 
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fg S r" log^ rAphg- ( |A.3| ) follows similarly from ( A.2 ) by Taylor expanding F. ( |A.4| ) and 



(A. 5) are elementary. ( A.6D follows by linearity from 



J r* log-' r dr 



log-'+-'^ r 
J + 1 



for some coefficients C, 



Appendix B Geodesies in polyhomogeneous metrics 

Proposition B.l ; Let g he a polyhomogeneous metric (Lorentzian or Riemannian) on a 
manifold with boundary M , g G Aphg n C^{M). 

1. For any p £ dM , k G TpM there exists e > and a unique geodesic Tp{s), s G [0, e), 
satisfying 

Tp{0)=p, tp{0) = k. (B.l) 
If we write Tp = {y'^{s)}, then y^{s) are polyhomogeneous functions of s. 

2. If dM B p ^ kp is a smooth field, k G C°^{dM), then the functions y^{s,v) are 
polyhomogeneous. 

Remark: For a polyhomogeneous metric we have dg ~ log^^ x near dM, so that standard 
results about geodesies do not apply. 

Proof: Let us define 

r{s) = y''{s)-sk>^-yi;, 

we thus have 

^ = F'^(V',V,s), (B.2) 

with 

F>^{i;, X, s) = r^^(yo" + sk-" + V) ■ (B.3) 

Let a G (0, |) and let e > be a number to be determined later. Consider the space 

X, = {^'',x''GC([0,e])} 

with the norm 

\\{^,x)\\x, = sup \s-'^x^{s)\+ sup Is— 

se[0,e] se[0,e] 

Let T : ^ X^he the map 

X, B = (^J\^{s)ds,Jy^{^P,x,s)ds^ , (B.4) 



F^^ given by (p. 31). Clearly a fixed point of T satisfies (B.l)-(p?3). From the estimates 



|r| + |5^,r| < ciog^x, (B.5) 

\d^r\ < Cx~^log^x, (B.6) 
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for some A^, it is easily seen that one can choose a constant K and an e > such that T 
is a contraction mapping from a bah around (0, 0) in of radius K into itself, and the 
result follows by the contraction mapping principle. Once the solution is known to exist, 
poly homogeneity immediately follows from the equation 

If A;^(f) is a smooth function of f G dM, then by considering the equations satisfied by 
d'^y^ (note that for a poly homogeneous metric d'^T satisfies the same estimates (|B.5| )~ ([B.6| ) 
as r itself) one obtains the result by an argument similar to the one above. □ 



Let us finally point out that part 1 of Proposition |B.l| is still true under the rather 
weaker hypotheses 



+ 15^,1 + x^-^\d^g^,\ + x^-^\d^dpg,,,\ < C , (B.7) 



for some constant C, with any /? > 0. The same proof as above goes through, except that 
the exponent a in the norm *° chosen to lie in (0, Gauss coordinates 



can be constructed for metrics satisfying ( B.7 ) provided that one moreover has 



x^-^\dAdpg^u\ + x^-^\dAdpd„g^,A < C , (B.8) 

where Ba are derivatives in directions tangent to dM, while 5p, etc., denotes all partial 
derivatives. 



Appendix C The Einstein field equations for the metric 
(2.1) 



Using SHEEP we have derived the Einstein equations for a metric of the form ( p. 11) . The 



metric coefficients have been further parametrized as in Proposition 3.2 with the small 
changes that, to aid comparison with is written as U and is written as W cosec 9. 

Derivatives with respect to the coordinates, which are numbered by 

{xq, xi, X2, X3) = {u, r, e, (p) 

are indicated by subscripts. The equations one obtains coincide with those of the Ap- 
pendix in [Q, except for some misprints listed below. If all terms were moved to the left 
hand sides, the left sides would be rRu/A in (C.l), 2r^Ri2 in ( |C.2 ), 2r^i?i3 cosec in 
0, -e^f^{N'^Rab)/2 in (Q, e^^{e-^^R22 - e^^ sm-^ 9 R^^) / Ar in O and e^'^{R2z - 
{h''^Rab)h23/2)/{2rcosli26sm9) in ([Oel) . 

The misprints in the corresponding equations of are: 

1. on p. 121, at the end of the second line of the second equation, a right parenthesis is 
missing; 

2. on p. 122, in the fourth line of the equation the signs of the terms which contain (323 
and (3203 should be reversed, and 



3. in the fifth line of the same equation the factor r should be replaced by 
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/3i = ^(7fcosh2 2,5 + ,5?) (C.l) 
,-2/3^g27^^ cosh 25 + Wi sinh 25)) ^ 
2r^(A2 + 2(5i(52 - 2r~V2 - ^1\52 cosh 2(5 sinh 2(5 
- (712 - 27172 + 271 cot 6) cosh^ 25) 

+2r'^e^^ cosec^^ — (5i3 — 2(5i73 + (713 + 27173) cosh 2(5 sinh 2(5 

+271(53(1 + 2 sinh^ 2(5)) (C.2) 



sinh 2(5 + 6-2^1^1 cosh 2(5) j ^ 
2r2e"^'^ ( - ^12 + 2(5i72 - 25i cot 6* - (712 - 27172 + 271 cot 9) cosh 25 sinh 25 
-271(52(1 + 2 sinh^ 2(5)) + 2r2 cosec 6 (/3i3 + 2(5i(53 - 2r-^l3^ 

+471(53 cosh 2(5 sinh 25 + (713 + 27173) cosh^ 2(5) (C.3) 

26^^^ cosec 6'((/323 + 1^20^ + 2(52(53) sinh 25 + ((523 + (^3 cot 6 + (5273 - i2h 
+<52/33 + /32<53) cosh 2(5) - e^'^-^T ((/322 + /3| + /32 cot e + 2-il + 2(5| - 1 

-722 - 372 cot e - 2/3272) cosh 2(5 + ((522 + 3(52 cot 6^ + 2/32-52 - 472-52) sinh 2(5) 
-6^/^+2^ cosec^ ((/?33 + /3| + 27I + 2<5i + 733 + 2/^373) cosh 2(5 + (-533 + 2^^ 

+473(53) sinh 2(5) - ((e^^C/i + e'^^'^W^) cosh 2(5 + 2[/iWi sinh 2(5) 

+^r (rC/12 + rC/i cot e + W2 + 4U cot ^) + ^r^ cosec ^ (VF13 + 4VF3) (C.4) 

01 cosh 2(5 + 2r(7o(5i + (5o7i) sinh 25 

UliVi + 711^ + r-^jiV) cosh 2(5 + 27i(5iy sinh 2(5 + Ir^e'^^ {e^^'Uf - e'^^W^) 
2 8 

+lr-'e'^-'^{(322 + ^1 - /32 cot 6) - ^^e'^+^^iPs^ + /3|) cosec^ 



+r-^e2'^ (/32(53 - P352) cosec + ^re'^'^ cosec 6 (^{Un + 2r-^C/3) sinh 2(5 

+A5iU3 cosh 2(5) - ^re'^T ( (H^i2 - Wi cot 0) sinh 2^ + 2r-^ {W2 - W cot e) sinh 25 

+4(5i (W^2 - W cot ^) cosh 2(5) - ([/12 + 2r-^?72 - Ui cot - 2r-^U col 

+4r~^'j2U + 47i2f/ + 272?7i + 271^72 + 271 [/ cot o) cosh 2(5 



-r {5iU2 + 271(52 [/ + 2(5i72f/ - SiUcot 9) sinh 2(5 

+-r cosec 6l(W"i3 + 2r"^H^3 - 4r"S3^ - 4713 - 273^^1 - 271^^3) cosh 2(5 
+r cosec 9 {5iWz - 25i'yzW - 2-fi5zW) sinh 2(5 (C.5) 



24 



{rd)oi — 2r7o7i sinh 26 cosh 26 
= ^{6iVi + 6nV + r-^6iV - 2jfV cosh 26 sinh 26) 

_lr-ie2/3~27(^22 + /3| - /?2 cot 6) sinh 2(5 - ir-ie2^+27 cosec^ 0(/333 + /^l) sinh 25 
-r-^e^'^ cosec (-/?23 - + Ps cot + /?273 - 72/33) cosh 26 

+lr^e-'^'^ Ue^^Ul + e-^^VF^^) sinh 25 + 2UiWi cosh2j) 
8 ^ ^ 

--r(2(5i2?7 + 2r-i52^7 + 61U2 + 52[/i + 5i?7cot 6 - 2{jiU2 - 71 C/ cot 6 + 47172^7) 

X cosh 25 sinh 25) - cosec ^ (2513VF + 2r-^53VF + 61W3 + 63W1 
+2 (71 VF3 - 27173 VF) cosh 25 sinh 25) - ^re'^'^ (Wu - Wi cot 9 
+2r~i {W2 - W cot 6) - 471 {W2 - W cot 6) cosh^ 25) - \re^^ cosec 9 



4 

x(C/i3 + 2r"^C/3 + 47iC/3Cosh2 25) (C.6) 

Appendix D Expansion in the axisymmetric case 

The foUowing formulae give the expansion of the axisymmetric case discussed at the end of 
section 3. Ah coefficients, except the last one in each quantity (i.e. 74, /^s, and V3) are 
considered to be functions of u and 9 only: 74, /?5, C/5 and V3 are written as functions of u, 
9 and r so that the consistency of the approximation can be checked by looking at the first 
neglected terms in the equations of appendix The exponentials of /3 and 7 are expanded 
to order r~^. 

As a result of the arguments of section 2, we know we can expand 7 as 

7 = c{v)r~^ + 72(w)r~^ + J3,i{v)r~^ logr + 73(v)r~^ + 74(r, f)r~'^ , 

where we use (v) to denote (n, x'*), and the dependence of 74 on r allows for logr terms 
there. To avoid confusion of subscripts in what follows we use the form f^x for the partial 
derivatives, where x is a variable name, not a number. 



Substitution in (|C.lD yields 



^ ,2^-2 2 3 3 ^-4,^„^ I ( ^ ^„,2\^-4 I n -5 



l3 = --c r - -C72r - -C73,ir logr + {^^^3,1 - |C73 - ^12 ] r + P^r 



Putting this into (C.2) gives 



U = -(2ccot(9 + ce)r"2- -72Cot6' + -72,0 r"^logr + r"3[73 



4^ 
3' 

+ (^2c72,e + 4c72 cot 9 + ^73,1,0 + 373,1 cot 9^ r"^ log r 

/5 5 5 3 1 2 

+ f -c^ cot 9 + -c^c e + -C72 cot 9 - -cUs + -072,0 + 272C ( 

11 11 3 \ _4 _5 

+y73,icot6' + 373Cot6'+ Y73,i,e + -73,ej r +r C/5 
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where U3 is as yet arbitrary. Its u-derivative is given by the -R02 Einstein equation, and on 
substituting van der Burg's form for the coefficient 

Us = 4c^ cot d + 3cce + 2iV 

we find agreement with his equation for A^^„. 
Next we use (p.4|) to obtain 



V = r - 2M + ^2 cot 6*72,0 - ^72 + ^72,e0^ r ^ log r 

+ (72,9 cot e + -^{cef + ^12,99 - 2^^3,5 + 4c^ cot^ e - -c^ - -72 
+-CC + ycc e cot - -C/3 cot 6* 1 r 

+ (8c^2 cot^ 6 + 4c72,6) cot 9 + 4720,9 cot 6* + 2c,6l72,6» 

13 \ 

- 2 73,1,061 - 2 73,1,0 cot 6* + 73,1 j r"^ log r 

/ o 2., 1 s 9 2 1 2 11 1 

+ 4c cot 6* + -c'' + -c c,6i cot fci - -c c^e9 + -77^2,9 cot fc* + -072,610 
V 4 4 b 

7 2 7 4 1 3 

+ o72C,0 cot 6 + -^20 99 + 7:0,912,9 + 7012 cot^ 9 - -C72 - ScUs cot 6 - -U3C,e 
6 6066 2 

5 15 ^51 3 . \ -2 

-g 73,1,00 - y 73,1,0 cot + -73,1 - -73,00 - -73,0 cot + 73 I r 

+r-^V3 

where M is as yet undetermined, but will have a u-derivative given by the i?oo equation as 

3 2 1 

M,u = - cot 9c,u9 - C,u - Cu + -C,u99- 

On substituting all these expansions into the equation ( p.5| ) we find that c,u is unde- 
termined, and 72,u = as expected. 73,i,n is then u-independent so we can integrate with 
respect to u and get 

73,1 = 73,1 + ^72,00 + ^72,0 cot 9 - ^72 cot^ 9 - ^72^ u 

where 73 is a freely specifiable function of 9. 

The next term in ( |C.5| ) gives the u-derivative of 73 as 

73,n = ^cc,e9 + ^cc,g cot 9 + ^{c,ef -c^cot'^ 9- l;c^ + l-cM 
00 o Z Z 

-^■72,99 - 77772,0 cot 9 + ^72 cot^ 9 - -^72 - ^U3,e + ^t/3 cot 9 

iZ iZ o 000 

and we cannot integrate this explicitly since we do not know the u dependence of c and 
hence of M and C/3. We can check, using 

to get van der Burg's form of the coefficient, that our result agrees with his when 72 = 0. 
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Finally we can solve the next order in (|C.5| ) for 74^^ and find that the result contains 
logr terms but no higher powers of logr, and that 

74,„ sin^ e = G + Fg 

where G is as in ( |3.12| ) and 



F = { ( -^72,6»06isin 61 + —72,00 cos 6* sin 6^ + —72,61 cos 6'sin( 

+ o72,9sin 0- -72 cos - -72Cos6'sin 1 
o 6 2 

--C72 cos 9 sin^ 6 — -072,6* sin'^ 9 — t(73,i),6» sin^ 9 + -73,1 cos 9 sin^ 9 | log r 

O O 4t Zi j 

+ \ -^12,996 sin^6' + ^12,89 COS 61 sin^ 6^ + ^72,0 cos^ 6lsin6' 

29 o 1 o 19 2 \ 

+^72,e sin 9 - -72 cos 6* - —72 cos 9 sin 9 j u 

lol q2 q1 q5 n 

— -73 61 sin 9 H — 73 cos 9 sin H — 72c sin H C72 sin 9 H — C72 cos 9 sin 

4 ' 2 3 ' 12 ' 6 

11 7 13 3 

-— cos 61 sin^ 9 - -c^c^ sin^ ^ + -cC/3 sin^ 6* - — (7° i),0 sin^ 9 + -7° ^ cos 9 sin^ 6* 

12 8 ' 4 16 ' ' 8 ' 

Since G contains no logr terms, it follows that Q as defined in ( |3.13| ) is conserved. (We 
may again note that if 72 = = 73,1 this agrees with [^].) 
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